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Abstract 

The standard setting of quantum computation for continuous problems uses deter- 
ministic queries and the only source of randomness for quantum algorithms is through 
measurement. Without loss of generality we may consider quantum algorithms which 



use only one measurement. This setting is related to the worst case setting on a classical 
computer in the sense that the number of qubits needed to solve a continuous problem 
must be at least equal to the logarithm of the worst case information complexity of this 
problem. Since the number of qubits must be finite, we cannot solve continuous prob- 
qh lems on a quantum computer with infinite worst case information complexity. This can 

even happen for continuous problems with small randomized complexity on a classical 
computer. A simple example is integration of bounded continuous functions. 

To overcome this bad property that limits the power of quantum computation for 
continuous problems, we study the quantum setting in which randomized queries are 
allowed. This type of query is used in Shor's algorithm. The quantum setting with 
randomized queries is related to the randomized classical setting in the sense that the 
number of qubits needed to solve a continuous problem must be at least equal to the 
logarithm of the randomized information complexity of this problem. Hence, there 
is also a limit to the power of the quantum setting with randomized queries since we 
cannot solve continuous problems with infinite randomized information complexity. An 
example is approximation of bounded continuous functions. 

We study the quantum setting with randomized queries for a number of prob- 
lems in terms of the query and qubit complexities defined as the minimal number 
of queries/qubits needed to solve the problem to within e by a quantum algorithm. 
We prove that for path integration we have an exponential improvement for the qubit 
complexity over the quantum setting with deterministic queries. 
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and the National Science Foundation. 
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1 Introduction 



One of the challenging problems of computational theory is the study of the power of quantum 
computation. By now there seems to be agreement about the standard setting of quantum 
computation, see [2, 7, 20]. This setting describes quantum computation as a sequence 
of unitary 2 k x 2 fc matrices acting on an initial state followed by a measurement. Here, k 
denotes the number of qubits. The unitary matrices in quantum computation are represented 
by elementary quantum gates, and one of them may represent a query which depends on the 
problem we want to solve. For continuous problems, queries are deterministic and depend 
on function values. Quantum algorithms may have many measurements but it is known 
that without loss of generality it is enough to consider quantum algorithms with only one 
measurement, see Remarks 1 and 2 and papers cited there. 

In what follows it is important to stress the difference between the cost of an algorithm 
for solving a given problem and the computational complexity of this problem. The com- 
putational complexity (for brevity, the complexity) is the minimal computational resources 
needed to solve the problem. Examples of computational resources which have been studied 
include memory, time, and communications on a classical computer, and qubits, quantum 
gates and queries on a quantum computer. 

We study quantum computation for continuous problems which are usually defined on 
spaces of functions. The quantum complexity of continuous problems has been studied in 
many papers and different queries, such as bit, phase and power queries, have been analyzed 
in the literature, see e.g., [4, 7, 8, 9, 10, 14, 15, 20, 21, 24, 26, 35]. In this paper we study bit 
queries, although some results hold for more general queries. What is important for our study 
is that a query for a continuous problem depends on at most 2 k function values computed at 
a priori given deterministic sample points. This means that in the standard quantum setting 
for continuous problems, quantum algorithms depend on at most 2 k function values and the 
only source of randomness comes through measurement. 

The query complexity has been the focus of research. It is defined as the minimal number 
of queries needed to solve a given problem to within e by a quantum algorithm. Since a 
critical resource for the foreseeable future is the number of qubits, we also study the qubit 
complexity which is defined as the minimal number of qubits needed to solve a given problem 
to within £ by a quantum algorithm. 

We stress that there could be a trade-off between the query and qubit complexities since 
the minimization of queries may lead to a large number of qubits and vice versa. We do 
not know of such trade-offs for continuous problems studied so far. It is unknown if such 
trade-offs occur for some continuous problems. 

To compute the quantum speedup one needs to know the worst case and randomized 
complexities on a classical computer. For continuous problems, the worst case and random- 
ized classical complexities have been thoroughly studied in information-based complexity, see 
[7, 21, 28, 32, 33, 38]. For our purpose, we need the concept of (non-adaptive) information 
complexity which is defined as the minimal number of function values needed to solve the 
problem to within e. We included two short sections on these classical settings to the extent 
needed in the rest of the paper. 

Our first technical result is a relation between the standard quantum setting and the worst 
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case classical setting. Namely, it is relatively easy to see that since quantum algorithms are 
based on at most 2 k function values, they can not have a quantum error smaller than the 
worst case error of a classical algorithm based on these 2 k function values. This analogy 
is not complete since in the worst case setting we use deterministic algorithms whereas a 
quantum algorithm has a random element through measurement. Nevertheless, it is possible 
to show that the qubit complexity is bounded from below by the logarithm 1 of the worst 
case information complexity of the problem which we want to solve to within 2e. We will 
show that this extra factor 2 takes care of randomness of quantum algorithms 2 . Since the 
worst case information complexity usually goes to infinity as e tends to zero, the number of 
qubits must also increase to infinity although at a much slower rate due to the presence of 
the logarithm. 

We also show that the qubit complexity is bounded from below by the (Kolmogorov) 
e-entropy of the solution set. Hence, problems with large entropy of the solution set require 
a large number of qubits. 

When the worst case information complexity or the entropy of the solution set is infinite 
then a finite number of qubits is not enough and the problem is unsolvable in the standard 
quantum setting. This can even happen for problems for which the randomized classical 
complexity is small. An example of such a problem is multivariate integration of continuous 
ci-variate functions defined on, say, [0, l] d , whose absolute values are bounded by 1. It is 
known that in this case the worst case information complexity is infinite but the Monte 
Carlo is optimal and the randomized information (as well as the total) complexity is roughly 
s~ 2 independent of d. 

Why can we solve this problem in the classical randomized setting and not in the standard 
quantum setting? The reason is that in the randomized setting we use function values at 
randomized points and potentially we can compute the function value at any point whereas 
in the standard quantum setting we use function values at deterministic points. The number 
of these points can be enormous, up to 2 k . But if we take a continuous function which 
vanishes at these 2 k points then we are unable to detect whether this function is zero or 
perhaps takes values equal to 1 or —1 at all points except an arbitrarily small neighborhoods 
of points at which it vanishes. The true solution may be zero or arbitrarily close to 1 or 
to —1. That is why any quantum algorithm in the standard quantum setting must fail. 

This negative result is our point of departure. To overcome this bad property of quan- 
tum algorithms and to enlarge the power of quantum computation we propose a small mod- 
ification of the standard quantum setting by allowing the use of randomized queries and 
randomized unitary matrices. The other assumptions are kept intact. We will call this modi- 
fication as the quantum setting with randomized queries and refer sometimes to the standard 
quantum setting as the quantum setting with deterministic queries. 

In fact, the idea of using randomized queries is not new. A particular kind of randomized 
query is used in Shor's algorithm for factoring of a (large) integer N, see [30] and also [20]. 
The essential part of Shor's algorithm is order finding which is solved by the query 

Qx\j) = \jx mod AO 

1 All logarithms in this paper are base 2. 

2 As indicated in the proof of Theorem 3.1 the extra factor 2 can be often omitted. 
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for j — 0,1, ... , 2^ log ^ — 1} with a random x from {2, 3, . . . , N — 1}. 

The use of randomized queries for continuous problems was also suggested in [24] for 
integration of non-smooth functions. 

The quantum setting with randomized queries is the same as the standard quantum 
setting with the one important exception that queries as well as all unitary matrices used 
by a quantum algorithm may now depend on a random element. Hence, we have now two 
sources of randomness: one affecting unitary matrices and the other affecting measurement. 

The quantum setting with randomized queries is an extension of the standard quantum 
setting. Indeed, if one always selects the same unitary matrices including the query, then 
we have exactly the standard quantum setting. Obviously, the use of randomized unitary 
matrices and randomized queries offers a possibility of much more efficient quantum com- 
putation. For some continuous problems, this extension is necessary. For example, we will 
show that multivariate integration of bounded and continuous functions, which cannot be 
solved in the standard quantum setting, is solvable in the quantum setting with randomized 
queries by a quantum algorithm that uses of order e _1 queries and log e^ 1 qubits. Hence, 
we have a quadratic speedup over the randomized setting on a classical computer. 

We now comment on the error criteria used in the quantum settings with deterministic 
and randomized queries. In the standard quantum setting, i.e., in the quantum setting with 
deterministic queries, two error criteria are studied: 

• the first error criterion is defined by taking the average performance of a quantum 
algorithm with respect to measurements for a worst function from the given class, 

• the second error criterion is defined by taking the worst case performance of a quantum 
algorithm with respect to measurements on a set of measure 1 — 5 for a worst function 
from the given class. 

The same error criteria are used in the quantum setting with randomized queries. In this 
case, randomization is richer and we take the average performance or a set of measure 1 — 5 
with respect to "measurements, randomized queries and randomized unitary matrices" . The 
first error criterion is studied in the main body of the paper whereas the second one is studied 
in the appendix. 

We define the randomized (bit) query and qubit complexities analogously to the ran- 
domized setting on a classical computer. The randomized query complexity is defined as 
the minimum of the average number of randomized queries needed to solve the problem to 
within e by a quantum algorithm. By "to within e" , we now mean that the error of a quan- 
tum algorithm is at most e which is defined by taking the average performance with respect 
to all random elements of the quantum algorithm for a worst function from the given class. 

The randomized qubit complexity is defined analogously as the minimal number of qubits 
for which there is a quantum algorithm whose error is at most e. We stress that we assume the 
number of qubits is fixed and does not vary during quantum computation. This is probably 
a reasonable assumption from a practical point of view since a quantum computer with a 
random number of qubits seems too much to be expected in the near future. Nevertheless, 
from a purely theoretical point of view it would be interesting to study also the quantum 
setting with randomized queries and with random number of qubits and try to minimize the 
average number of qubits needed to solve the problem. 
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It is not surprising that the quantum setting with randomized queries is related to the 
randomized setting on a classical computer in the sense that the randomized qubit complex- 
ity is bounded from below by the logarithm of the randomized (non- adaptive) information 
complexity on a classical computer. If the randomized information complexity of a problem 
is infinite, the problem cannot be solved in the quantum setting with randomized queries. 
This happens, for example, for approximation of bounded continuous functions. Of course, 
the class of problems with infinite randomized information complexity is smaller that the 
class of problems with infinite worst case information complexity. So we extend the limit of 
what can be computed by presenting the quantum setting with randomized queries. 

We study the quantum setting with randomized queries for a number of problems, and 
for some of them we prove an exponential improvement for the qubit complexity compared 
to the standard quantum setting. This is especially important since, as already mentioned, 
the number of qubits is a critical resource for the foreseeable future. In particular, the 
exponential improvement holds for path integration. 

In this paper, we study real and Boolean summation, multivariate integration and path 
integration. We now briefly state the results obtained for these problems. 

The real summation problem lies at the core of many continuous algorithms and plays a 
major role in the study of continuous problems in the standard quantum setting. The same 
is true in the quantum setting with randomized queries. It is known, see e.g., [7, 21], that the 
real summation problem can be reduced to Boolean summation. That is why it is enough 
to present in detail results for only the latter problem in which we want to approximate 

JV-l 
3=0 

for a Boolean function / : {0, 1, . . . , N — 1} — > {0, 1}. Here N is a large integer which can 
be assumed to be a power of 2. We want to compute B^if) to within e. Without loss of 
generality we may consider s~ l <C N. We now present the orders of the query and qubit 
complexities in the quantum settings with deterministic and randomized queries. 





Deterministic Queries 


Randomized Queries 


Query Complexity 


s- 1 


e- 1 


Qubit Complexity 


log N 


log e^ 1 



Figure 1: Boolean Summation 



We stress that the minimal numbers of queries and qubits in a given setting are obtained by 
essentially the same quantum algorithm. In the quantum setting with deterministic queries, 
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this is the Boolean summation algorithm of [5] with seven repetitions as proved in [11]. In 
the quantum setting with randomized queries, we first approximate the Boolean mean Bjy(f) 
by the Monte Carlo algorithm, 

1 m 

m A — ' 

with m of order e~ 2 and with independent and uniformly distributed u>j from {0, 1, . . . , N — 
1}, and then use the Boolean summation algorithm with seven repetitions to approximate 
MC m (/,w). It is interesting to notice that this algorithm uses randomized queries but the 
remaining unitary matrices are deterministic. This leads to the corresponding upper bounds. 
The lower bound proof of the query complexity in the standard case follows from [18]. For 
the quantum setting with randomized queries, we use the known fact that the randomized 
errors of quantum algorithms are no smaller than the average case errors with respect to 
Boolean functions. The latter problem with an appropriate measure on Boolean functions 
was solved in [25]. The lower bound proof of the qubit complexity is from the relation to 
the randomized information complexity. 

We stress that we have the same order of query complexities in both cases. However, by 
allowing randomized queries we have an essential improvement in the number of qubits for 
solving the Boolean summation problem. 

For the real summation problem, we want to approximate 

N-l 

SUM„(/) = - 

3=0 

where / : {0, 1, . . . , TV — 1} — > [0, 1] may now take real values. For completeness, in Section 5 
we show how the real summation problem may be reduced to the Boolean summation prob- 
lem. In Corollary 5.1 we show that the results presented in Figure 1 also hold for the real 
summation problem. 

We now turn to multivariate integration for functions which are r times differentiable and 
uniformly bounded. For r = 0, the query and qubit complexities are infinity in the quantum 
setting with deterministic queries. For the quantum setting with randomized queries, they 
are finite and their orders are given in the following table. 





Deterministic Queries 


Randomized Queries 


Query Complexity 


oo 


e- 1 


Qubit Complexity 


oo 


log e~ l 



Figure 2: Multivariate Integration for r = 
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Hence, in this case the improvement of the quantum setting with randomized queries 
over the standard quantum setting is infinite. 

We now assume that r > 1. Hence, functions are now at least once differentiable. In this 
case the orders of the query and qubits complexities are the same in both cases. 





Deterministic Queries 


Randomized Queries 


Query Complexity 


£ -l/(l+r/d) 


e -l/(l+r/d) 


Qubit Complexity 


log E^ 1 


log e^ 1 



Figure 3: Multivariate Integration for r > 1 



The query complexity in the standard quantum setting is due to [24]. The randomized 
query complexity has the same order since Boolean and real summation require roughly 
the same queries in both settings. The qubit complexities are of the same order since 
the logarithms of the worst case and randomized information complexities of multivariate 
integration are both proportional to log e^ 1 . 

Finally, we consider a specific case of path integration studied in [35]. The orders of 
query and qubit complexities are presented in the following table. 





Deterministic Queries 


Randomized Queries 


Query Complexity 


e -i+°(i) 


e -i+°(i) 


Qubit Complexity 


s~ 2 log e^ 1 


log e^ 1 



Figure 4: Path Integration 



We thus have the same orders of query complexities and an exponential improvement in 
the number of qubits. 

We stress that in the randomized classical setting and in the quantum setting with 
randomized queries we permit the use of random elements from a set Q whose cardinality 
may be infinite and distribution of points from Q may be arbitrary. For example, the 
classical Monte Carlo with n random points for integration of c?-variate functions defined 
over, say, [0, l] d , uses Q = [0,1] dn and uniform distribution. Alternatively, it is possible, 
also for classical computers, to use a restricted form of randomization based on, for example, 
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random bits or a finite set f2. This obviously restricts the class of randomized algorithms and 
it is not clear if positive results for unrestrictive randomization are still true for the restricted 
case. There is a very interesting stream of work, see [13, 27], studying the minimal number 
of random bits needed for the solution of continuous problems on a classical computer. 
There are also general results in [12] showing that as long as Q is finite then the classical 
randomized setting is (roughly) equivalent to the standard quantum setting at the expense of 
adding additional qubits. However, if the cardinality of Q goes to infinity then the additional 
number of qubits also goes to infinity. That is why, the standard quantum setting is not 
equivalent to the classical randomized setting without a restriction on Q. 

We hope that the quantum setting with (restricted) randomized queries will be studied 
for general continuous problems. It would be especially interesting to characterize continuous 
problems for which this setting offers an exponential improvement in the number of queries 
and/or qubits over the standard quantum setting. 



2 Continuous Problems 

The computational complexity of approximate solutions of continuous problems has been 
studied in information-based complexity, see e.g., [7, 21, 28, 32, 33, 38]. We present a brief 
outline of this theory in the worst case, randomized and quantum settings to the extent 
needed for this paper. 

Let F be a non-empty subset of a linear space of o?-variate functions / : — > R with 
Dd C M. d . Let G be a normed space with its norm denoted by || • ||. Consider a (linear or 
non-linear) operator 

S : F -> G. 
Our goal is to compute S(f) to within e for / £ F. 

Example : Multivariate Integration 

We illustrate the concepts of this paper by an example of multivariate integration of smooth 
functions. Let C r ([0,l] d ) denote the class of real functions defined on the (/-dimensional 
unit cube, / : [0, l] d — > R, all of whose partial derivatives up to order r exist and are 
continuous. That is, for a multi-index a — [a±, a 2 , ■ ■ ■ , otd] with non-negative integers otj and 
with \a\ := ai + a 2 + ■ ■ ■ + ad < r we know that 

D a f = f 

J dt? 1 dt? 2 ---dt% d J 

exists and is continuous. The norm in C r ([0, l] d ) is defined as 

||/|| r = max max |£>"/(f)|. 

a: \a\<r te[0,l] d 

Then we set 

F = F d , r = {fEC r ([0,l} d ): \\f\\ r <l} 
as the unit ball of C r ([0, l] d ), and G = R. 
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The multivariate integration problem S = INT^r : C r ([0, l] d ) — > R is defined by 

INT d>r (/) = / /(t)cft. 

J[0,l] d 

This is an example of a linear problem since the operator INT djr depends linearly on /. □ 

The approximate computation of S(f) can be done as follows. First of all, we specify how 
information about the function / is used by algorithms. We assume that we can compute 
finitely many function values 3 f(t) for some sample points t from That is, any algorithm 
may use f(ti), f(t 2 ), ■ ■ ■ , f(t n ) for some n and tj. We stress that the choice of the sample 
points tj may be adaptive, i.e., t\ is given a priori, whereas tj may depend on the already 
computed values f(ti), ffe), ■ ■ ■ , f(tj-i) f° r j = 2, 3, . . . , n. The number n can also be chosen 
adaptively. The sample points may be chosen deterministically or randomly depending on 
the setting, for details see [32]. The specific form of an algorithm also depends on the setting 
in which we define the error. We first present two settings for classical computers and then 
turn our attention to the quantum setting with deterministic and randomized queries. 

2.1 Classical Computers: Worst Case Setting 

In the worst case setting, we assume that sample points as well as algorithms are determin- 
istic. That is, an algorithm that uses n function values has the form 

Mf) = 0(/(ii),/(t 2 ),...,/(y) (i) 

for some mapping : IR n — > G. If the sample points are given a priori and are the same for all 
/ from F, then A n uses non-adaptive information. Otherwise, it uses adaptive information. 

The worst case error of the algorithm A n is given by its worst case performance with 
respect to /, 

e™{A n ) = sup\\S(f)-A n (f)\\. 

feF 

Example : Multivariate Integration (continued) 

A typical choice of an algorithm for multivariate integration is a linear algorithm, sometimes 
called a quadrature or cubature, 

n 

Mf) = E%M) 

for some aj G R and tj G [0, l] d . 

For r > 1, it was proven by Bakhvalov already in 1959, see [1] and also [21, 32], that 
the minimal worst case error of algorithms using n function values is proportional to n~ r l d . 
Furthermore, the error of order n~ r / d can be achieved by a linear algorithm using non- 
adaptive information. Hence, if we want to guarantee that e WOT (A n ) < e, then n has to be of 

3 More general information given by arbitrary linear Junctionals on / has also been extensively studied in 
information-based complexity, see e.g., [32] 
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order e~ d ^ r , and is exponential in d. This is called the curse of dimensionality meaning that 
multivariate integration is intractable in the worst case setting if d is much larger that r. 
The case r = will be considered later. 

For large d, a popular choice of cij is n _1 which leads to QMC (quasi-Monte Carlo) 
algorithms. The sample points are often chosen as low discrepancy points, lattice or shifted 
lattice points, see [19, 29]. For some spaces other than C r ([0, l] d ) the error behavior of such 
algorithms is only polynomial in d or even independent of d and tends to zero as a positive 
power of n" 1 . This is an active research area of information-based complexity dealing with 
high dimensional problems; the reader may consult [23] for a survey. □ 



2.2 Classical Computers: Randomized Setting 

In the randomized setting, we allow randomized choices of sample points as well as algo- 
rithms. That is, we have a probability space of elements uo from some set Q which are 
distributed according to some probability measure p on Q, p(fi) = 1. Algorithms using n 
function values on the average have now the form 

A n (f; u ) = ^ (f(t^), /(f 2>w ), . . . , f(t n ^)) , (2) 

where tj tU! are randomized sample points from Dj, and W is a randomized mapping from 
M. n " to G. Here, is the randomized number of sample points and its average is n, i.e., 




We stress that the sample points tj :U! as well as n w can be chosen adaptively as in the worst 
case setting. This also means that the probability measure p may depend on the function 
/ through its computed function values, see Chapter 10 of [32] for details. If the sample 
points tj jU are the same for all / from F, then A n uses non-adaptive randomized information, 
otherwise it uses adaptive randomized information. 

The randomized error of the algorithm A n is defined by its worst case performance with 
respect to / and the average performance with respect to cu, 

e™(A n ) = sup ( f \\S(f) - A n {f,u)fp{du) \ ' . (3) 
feF \Jn J 

Here, we choose to study the average performance in the L 2 -norm; however it is also possible 
to study it in a more general case of the L p -norms with p e [1, oo). 

Example : Multivariate Integration (continued) 

Probably the most popular and widely used randomized algorithm is the Monte Carlo algo- 
rithm 

A n (i» := MC n (i» = -J2f( tjiUJ ), 

H 3=1 

where tj jU) are independent and uniformly distributed sample points over [0, l] d . In this case, 
Q = [0, l] dn and p is Lebesgue's measure. That is, uj = [cui, 002, ■ ■ ■ , u n ] with Uj G [0, l] d 
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and tj :Ul = ujj. We stress that Monte Carlo uses non- adaptive randomized information with 
deterministic n and the deterministic mapping W = given by <f>(y) = rT x Y^j=iVi- ^ ls 
well known that 

J[0,1]^ \J{o,i]* n ~t J n 

Since for / e with r > 0, we have Jj ^ f 2 (t)dt < 1, then 

e ran (MC n ) < n- 1/2 . 

Hence, e ran (MC n ) < e for n = \e~ 2 ] and the curse of dimensionality of the worst case setting 
is broken by Monte Carlo in the randomized setting. Bakhvalov also proved, see [1], that 
the minimal randomized error of algorithms using n function values is of order n~ 1 / 2+r / d , 
and the latter error bound is achieved by a linear algorithm using non-adaptive information. 
This bound also holds if we use n function values on the average as proven by Novak in [22]. 
Thus, Monte Carlo almost minimizes the randomized error if d is much larger than r. □ 

The errors of randomized algorithms may be defined differently than (3). This corre- 
sponds to the probabilistic errors which are related to the quantum setting error commonly 
used in many papers. To simplify the presentation of the paper, we deal with the probabilistic 
errors in the appendix. 

2.3 Complexity and Information Complexity 

As already mentioned, we want to compute S(f) to within e. That is, we are looking for an 
algorithm A n whose error in the worst case or randomized setting is at most e, 

e wor/ran (A n ) < e. (4) 

We would like to guarantee (4) with the minimal cost of computing A n (f). This minimal 
cost is called the (total) e-complexity of S, and denoted by comp WOT ^ n (e,S). The cost of 
computing y = A n (f) is defined by counting the cost of n function values plus all operations 
needed to obtain y. The abstraction typically used in information-based complexity (and in 
scientific computation) is the real number model of computation in which we assume we can 
perform arithmetic operations and comparisons of real numbers with unit cost independently 
of the size of numbers, again see [32, 33] for details. The reader is referred to [31] for 
the motivation behind the real number model and comparison with the Turing model of 
computation. 

As we shall see, for quantum computation the complexity of S is less relevant than the 
non-adaptive information complexity. The latter is defined as the minimal number n of 
non-adaptive function values in a given setting needed to find an algorithm A° ad with error 
at most e. More precisely, in the worst case setting, algorithms A° ad are of the form (1) with 
a priori given sample points tj, whereas in the randomized setting, they are of the form (2) 
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with sample points tj :U! independent of / and depending only on uj with fixed n w = n. Hence, 
the non-adaptive information complexity is defined by 

comp inf ~ wor / ran (e, S) = min { n : 3 Af d such that e wor / ran (^ ad ) < £ } . (5) 

We stress that the minimum in (5) is taken over all algorithms using non-adaptive infor- 
mation. That is, over all possible sample points tj and functions in (1) in the worst case 
setting, and over all probability measures p, sample points tj jU and functions W in (2) in the 
randomized setting. 

Surprisingly, for many continuous problems the non-adaptive information complexity is 
practically the same as the total complexity. There are, however, continuous problems for 
which the use of adaptive information is crucial and the non-adaptive information complexity 
is much larger than the total complexity, see [34] pp. 165-170. There are also continuous 
problems for which the reverse is true. That is, non-adaptive information complexity is small 
but it is impossible to combine it in a finite number of operations, and therefore the total 
complexity is infinite, see [36]. 

Example : Multivariate Integration (continued) 

Bakhvalov's results mean that the total and non-adaptive information complexities of mul- 
tivariate integration for the unit ball of C r ([0, l] d ) are of order 

e ~ d / r j n the worst case setting with r > 1, 
£ - 2 /( 1 + 2r / d ) j n the randomized setting with r > 0. 

□ 



3 Quantum Setting with Deterministic Queries 

We describe the standard quantum setting by presenting a general form of quantum algo- 
rithms used in this setting and the definition of their errors. Quantum algorithms can be 
characterized, in particular, by the number of queries and qubits they use. If they use n 
queries and k qubits, we will denote them by A n ^. Queries and qubits are deterministic and 
the only source of randomness is through measurement. We stress that quantum algorithms 
use non- adaptive information about the functions /. 

The quantum algorithms A n ^ are of the following form, see [5, 7, 21]. All computations 
are done on unit vectors in the complex space C . Here, k denotes the number of qubits. 
We assume that the initial state is a unit vector \ip ) from C 2 \ For / e F, the final state 
\i])f) is equal to 

= UnQfUn-iQf ••• £/iQ/ Wo), (6) 

where Uq, Ui, . . . , U n are 2 k x 2 k unitary matrices which are independent of /. Usually it is 
required that each Uj is represented by a relatively small number of elementary quantum 
gates. This will not be important for our considerations and we permit the use of arbitrary 
unitary matrices Uj. Of course, this makes lower bounds on the number of needed queries 
and qubits stronger. 
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The query Qf is also a 2 k x 2 k unitary matrix and depends on the function /. We assume 
in this paper that Qf is a bit query, see [7, 8, 9, 10, 20, 24] although results on the qubit 
complexity also hold for more general queries such as phase and power queries studied in 
[4,26]. 

For a Boolean function / : {0, 1, . . . , 2 m — 1} — > {0, 1}, with k — m + 1, the bit query is 
defined by 

Qf\m = \j)\i®fU)) 

for all i G {0, 1} and j G {0, 1, . . . , 2 m — 1}, and © denotes the addition modulo 2. 

For real functions / : D d — > R, we assume 4 that k = mi + m 2 , where mi qubits are 
needed to code the arguments of / and m 2 qubits are used for the values of /. The coding 
is done by two mappings 

r : {0,l,...,2 mi -l} - D d , 

P:f(D d ) -> {0,l,...,2 m2 -l}, 

and the bit query takes the form 

QAM) = \j)\i®P(f(T(j))) 

for all j G {0, 1, ... , 2 mi — 1} and % G {0, 1, ... , 2 m ' 2 — 1}, and © now means the addition 
modulo 2™ 2 , for more details see [7] . 

We stress that the bit query depends on at most 2 mi function values computed at some 
non-adaptive points tj = t(J). Furthermore, although we will not use this fact later, these 
function values are usually computed with some noise due to the finite range of the coding 
function (3. Usually, fl(f(tj)) is defined as the m 2 most significant bits of f(tj). Obviously, 
2 mi < 2 k . If f(D d ) is bounded then mi and m 2 are usually of the same order. To simplify 
further considerations we will use 2 k instead of 2 mi . 

For our purpose, the most important property of the bit query is that Qf depends on 
at most 2 k function values taken at some a priori given (non-adaptive) deterministic sample 
points tj from D d , 

Qf = Qf(t!),f(t 2 ),...,f(t 2k )- (7) 

The results on the qubit complexity will be derived using the property (7). Therefore 
they will be valid for all queries satisfying (7) which hold, in particular, for bit, phase and 
power queries. 

The bit query Qf is a deterministic 2 k x 2 k unitary matrix, and therefore the final state 
\ipf) is also a deterministic vector of 2 k components which use n times the query Qf based 
on non-adaptive information consisting of at most 2 k function values at some sample points. 
This means that if we consider two functions f\ and ji both from F such that fi(tj) = f2{tj) 
for j = 1, 2, . . . , 2 k then the queries Qf x and Qf 2 are the same, and therefore we obtain the 
same final states 

\M = \M 

for both fi and / 2 . 

4 For simplicity we do not consider ancilla qubits. 
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The only source of randomness is through measurement. That is, we obtain an index j e 
{0, 1, . . . , 2 k — 1} with probability pfj which depends on the final state \ipf). As before, pfj 



for some mapping : {0, 1, . . . , 2 fc — 1} — > G. The algorithm A ni fe is called a quantum 
algorithm. 

We stress that the quantum algorithm that uses fc qubits takes at most 2 k different values 
independently of the number n of queries used. 

Remark 1: We add that, in principle, we may use hybrid algorithms that are combinations 
of classical algorithms using bit operations on classical computers and quantum algorithms 
with many measurements. However, it is known, see e.g., [3, 7], that such algorithms may 
be rewritten in the form (6) and (8) with one measurement at the end of computation by 
linearly increasing the number of queries and qubits. It is also known, see Lemma 1 in [8], 
that we can sample r(/) instead of / if r(/) depends on k function values of /. Then the 
query Qr(f) can be simulated by a quantum algorithm that uses 2k queries on /. Therefore, 
without loss of generality we may consider only quantum algorithms with one measurement 
of the form (6) and (8). 

We stress that this is true if hybrid algorithms use only bit operations on classical com- 
puters. If we use the real number model of computation then not every algorithm can be 
written in the form (6) and (8). One reason of this is that we may have infinitely many out- 
puts in the real number model of computation which is impossible to obtain in the quantum 
setting. □ 

We now discuss the error of a quantum algorithm. There are at least two natural ways to 
define the error. One of them is by taking the worst case performance of a quantum algorithm 
with respect to / and the average case performance with respect to the index j. The other 
is to take the worst case performance with respect to / and the worst case performance with 
respect to the index j modulo a set of measure S for some (usually small) positive 5. For 
some problems, when S is a linear functional, it is enough to take, say, 5 = 3/4 and increase 
the probability of success by running the quantum algorithm a couple of times and by taking 
the median as the final result. 

We will study both definitions of the error of a quantum algorithm. In the main body of 
the paper we choose the first option since it is directly related to the error usually studied 
in the randomized classical setting. The other error, which is probably more popular in the 
quantum literature, is called the probabilistic error and is studied in the appendix. 

Hence, by the error of the quantum algorithm A n ^ we mean 



This concludes the definition of the standard quantum setting which can be summarized by 
the general form of a quantum algorithm (6) and its error (9). 




qua— std 



(A n>k ) = sup 
feF 
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We are interested in finding quantum algorithms with error at most e. We would like to 
achieve this goal with the minimal number of queries and/or qubits. 
By the query complexity in the standard quantum setting we mean 

comp que - std (£, S) = mm{n : 3A n , k such that e qua - std (A n , fc ) < e } . (10) 

By "there exists A n ^ we mean a quantum algorithm using n queries and k qubits with a 
finite k which can be, however, arbitrarily large. Hence, it may happen that the minimization 
of the number of queries will be possible at the expense of the number of qubits. 

By the qubit complexity in the standard quantum setting we mean 

comp qub - std ( £ , S) = min { k : 3^ such that e qua - std (A n ,fc) < e } . (11) 

In this case, by "there exists A n ^ we mean an arbitrary choice of the number n of queries. 
Obviously, k in (10) must be at least as large as comp qub ~ std (e, S), and n in (11) must be at 
least as large as comp quc ~ std (e, S). 

Although we do not pursue this point in the paper, it is also reasonable to minimize both 
queries and qubits. For example, we may want to minimize kn or n + [3 k, for a given positive 
number (3, over all quantum algorithms using n queries and k qubits whose quantum error 
is at most e. Here, if we choose (3 small then our emphasis will be on the number of queries, 
and if (3 is large then our emphasis will be on the number of qubits. 

Remark 2: We stress that query and qubit complexities are defined by minimizing the 
number of queries/qubits needed to solve the problem by a quantum algorithm with one 
measurement. 

Suppose we have quantum computation which requires the use of a sequence of quantum 
algorithms A n .^. each with one measurement and uses rij queries and kj qubits for j = 
1,2, ... ,p. Then the total number of queries is n = Yl P j=i n j which is of the same order when 
all A transformed as a quantum algorithm A with one measurement which uses 0{n) 

queries. 

The situation is, however, different for qubits since to run all quantum algorithms A nj ^, j 
it is enough to have k = maxj = i 5 2,..., p kj qubits whereas the quantum algorithm A would 
require of order Y^j=i kj qubits. Obviously, as long as p does not depend on e, it does not 
really matter since k must be at least of the same order as the qubit complexity. If, however, 
p is large and depends on e, then the qubit complexity needed for quantum algorithms with 
one measurement may be improved by many measurements. 

There is one case for which the size of p does not matter. Namely, when the qubit 
complexity is infinite which can happen as we see in the next section. 

It would be tempting to redefine the qubit complexity as the minimal number of qubits 
needed to solve the problem by a hybrid algorithm which performs classical and quantum 
operations with possible many measurements. This minimum is, however, zero since we could 
simulate all quantum operations on a classical computer with no qubits but at exponential 
cost of classical operations. 

We choose to study the qubit complexity of quantum algorithms with one measurement 
to eliminate such a possibility. □ 
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Example : Multivariate Integration (continued) 

For the class F = F d:T with r > 1, the minimal error e qua_std (A ri fc ) of quantum algorithms 
A ri:k is of order n~ 1 ~ r / d and is achieved by an algorithm that uses of order log e~ l qubits. 
This result follows from reduction of the integration problem to real and then to Boolean 
summation, and from the fact that the Boolean summation algorithm of [5] using n bit 
queries with seven repetitions has the error for worst / and average j also proportional to 
rT x as proved in [11]. This implies that 

comp que - std (£, S) = 6 (e-VCi+r/d)) and C omp qub - std (£, S) = O (log e' 1 ) . (12) 

□ 

3.1 Lower Bounds on Qubit Complexity 

We now prove lower bounds on the qubit complexity in the standard quantum setting in 
terms of the non-adaptive information complexity in the worst case setting as well as in 
terms of the (Kolmogorov) ^-entropy of the set S(F). Based on these bounds, we conclude 
that some continuous problems S cannot be solved in the standard quantum setting. 

Theorem 3.1. 

comp quh - std (e,S) > log comp inf - wor (2e, 1 S). 

Proof. 

Take an arbitrary quantum algorithm A n ^ such that e qua ~ std (A„ 5 fc) < e with the minimal 
number of qubits k = comp qub ~ std (£, S). We have 

The final state as well as probabilities pfj of the quantum algorithm A n ^ are based on 
the non-adaptive information 

N(f) = [/(fi), /(fa),---, /(**)] 
for some sample points tj G Therefore we can write 

An, k (fJ) = Hf(ti)J(t 2 ),...J(t 2k );j) VfEF, VjG {0,l,...,2 fc -l}, 

for some mapping $ : R 2fc x {0, 1, . . . , 2 k — 1} — > G. 

For an arbitrary / G F, take two functions fx and f 2 such that N(fi) = N(f 2 ) = N(f). 
The final state as well as all probabilities pfj will be the same for f\ and f 2 , and aj = 
A n ,k(fiJ) = Ai,k{f^j) for a11 3- Hence, for any / G F, we have 

2^ 2 > J^Pf,j (\\S(h) - %|| 2 + ||5(/ 2 ) - % || 2 ) . 

3=0 
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Observe that 

Wi)-S(/ 2 )|| 2 < (\\S(fi)-a i \\ + \\S(f 2 )-a j \\) 2 <2(\\S(f 1 )-a j \\ 2 + \\S(f 2 )-a j \\ 2 ). 
Multiplying both sides by pfj and summing up with respect to j, we conclude 

IWi)-S(/ 2 )H 2 <4e 2 . 
Taking the supremum with respect to / e F and /i, / 2 from F with N(fi) = N(f 2 ) we have 
sup sup ||5(A) - 5(/ 2 )|| < 2e. 

feF f 1 J 2 eF,N(f 1 )=N(f 2 )=N(f) 

The left-hand side of the last inequality is equal to the diameter of information N, see 
[32] p. 45, which in turn is bounded from below by the radius of information, denoted by 
rad(iV). Hence, rad(iV) < 2e which can hold only if the cardinality of N is at least equal to 
comp inf - wor (2e,S), see [ 32 ] p 54. ThuS; 2 k > comp inf - wor (2e, S), as claimed. 

We add in passing that for many cases we have d(N) = 2rad(iV). This holds, in particular, 
if G — R. Then rad(iV) < e and the extra factor 2 can be omitted. □ 

Theorem 3.1 states that the number of qubits needed to solve S in the standard quantum 
setting is related to the non-adaptive information complexity in the worst case setting. For 
most continuous problems S, the non-adaptive information complexity goes to infinity as 
e approaches zero. Then Theorem 3.1 says that the number of qubits also goes to infinity 
although much more slower due to the presence of the logarithm. We illustrate this point 
by continuing our example. 

Example : Multivariate Integration (continued) 

Consider F = F^ r with r > 1. We know that 

comp inf - wor (2e, 1 S) = 6 (e~ dlr ) . 
Then Theorem 3.1 supplies a lower bound on the qubit complexity, 

comp qub - std (£, INT d>r ) > f \oge~ 1 + 

with the term fl(l) independent of e but dependent on d and r. 

As we already discussed, - was the exponent of the worst case complexity of the integra- 
tion problem INT d r and caused the curse of dimensionality. Its role for the qubit complexity 
is mitigated since it effects the lower bound on the qubit complexity only linearly. 

Due to (12) the lower bound on the qubit complexity is sharp with respect to e, and we 
have 

comp qub - std (e,INT d , r ) = (log s' 1 ) . 

The dependence on e is very weak although for e tending to zero, the qubit complexity slowly 
goes to infinity. □ 

Remark 3: Theorem 3.1 was formally proved for the class of quantum algorithms with one 
measurement. We now show that a similar result holds for the much more larger class of 
hybrid algorithms which use non-adaptive or adaptive function values on a classical computer 
and many measurements on a quantum computer. More precisely, consider the following class 
of hybrid algorithms: 
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• For i — 1,2, ... ,p do 

— Use a classical algorithm with TOj non-adaptive or adaptive function values to get 
an initial state \ipo,i), 

— Use a quantum algorithm A n .^ with one measurement starting with the initial 
state iV'o.i) an d with bit queries and fcj qubits, and let 

A n . M (f,j) = <Pi{j) 

with a function 0j which can now be dependent on ^ non-adaptive or adaptive 
function values. 

Observe that the total number of function values used by the hybrid algorithm is at most 

p 

i=i 

and up to Y^=i( m i + ^) °f them can be computed adaptively. The hybrid algorithm uses k 
qubits, where 

k = max fcj. 

i=l,2,. ..,p 

Assume that the error of the hybrid algorithm is e. Then as in the proof of Theorem 3.1 we 
conclude 

N > comp inf - wor (2e,5), 

where now comp mf ~ wor (2£, S) stands for the worst case (adaptive) information complexity 
defined as in (5) with the exception that now A^ ad is replaced by an arbitrary algorithm A n 
using at most n adaptive function values. 

Hence, as long as there are two positive numbers ai and a 2 such that 

N < a 1 2 a2k 

then 

k > ajMog comp inf - wor (2e, 1 S) - log a x . 

Hence, even for hybrid algorithms, the logarithm of the worst case (adaptive) information 
complexity is essential and tells us how many qubits are needed. □ 

We now consider the case when comp inf_wor (£, S) = oo, i.e., when we cannot solve the 
problem in the worst case setting. Then the qubit complexity is also infinite. We summarize 
this fact in the following corollary. 

Corollary 3.1. // the non-adaptive information complexity of S in the worst case setting is 
infinity then S cannot be solved in the standard quantum setting. 

We illustrate Corollary 3.1 by multivariate integration for r = 0. 
Example : Multivariate Integration (continued) 

Assume now that r = 0. Hence, F = F^o is the unit ball of continuous functions with the 
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norm ||/||o = m &x- x e[o,i] d 1/(^)1 bounded by one. It is known, and easy to see, that for any 
algorithm A n we have 

e wor (A n ) > 1 Vn. 

Indeed, as already explained in the introduction, it is enough to take two continuous functions 
from F vanishing at the sample points tj used by the algorithm A n , j = 1,2, ... ,n, such 
that the integral of the first function is almost 1, and the integral of the other function is 
almost —1. Since these functions are indistinguishable for the algorithm A n , the best we can 
do is to approximate their integrals by zero with error arbitrarily close to one. Hence, the 
worst case error of any algorithm is at least one, as claimed. This implies that 

comp inf - wor (e,5) = oo Ve < 1. 

Theorem 3.1 says that multivariate integration for r = is unsolvable in the standard 
quantum setting. 

As already mentioned, this problem is, however, easy in the randomized setting. The ran- 
domized error of the Monte Carlo is bounded by rT 1 ! 2 which is optimal due to lower bounds 
of Bakhvalov and Novak. Therefore the non-adaptive randomized information complexity as 
well as the total randomized complexity are both of order e^ 1 ^ 2 . 

So why does the standard quantum setting fail for the problem which is relatively easy 
in the randomized setting? As we shall see in the next section, the reason is that we use 
deterministic queries in the standard quantum setting. This bad property will disappear if 
we allow the use of randomized queries also in the quantum setting. □ 

Before we proceed to the quantum setting with randomized queries, we briefly present 
another lower bound on the qubit complexity in the standard quantum setting. This bound 
relates the qubit complexity to the (Kolmogorov) e-entropy of the set S(F). We first recall 
the notion of e-entropy in normed spaces, see e.g., [17]. Let B be a subset of G. We want to 
cover the subset B by the minimal number of subsets of G whose diameters do not exceed 2s. 
That is, let 

n(e, B) = min { n : 3 Bj C G such that diam( J B i ) <2s, Be U] =1 Bj } , 
where diam(Sj) = sup bl b2eBj \\bi — 6 2 ||. Then the e-entropy of B is 

Ent(e,S) = log n(e,B). 
It is easy to prove the following theorem. 
Theorem 3.2. 

comp qub - std (£,S) > Ent(£,S(F)). 

Proof. 

The proof relies on the fact that in the standard quantum setting any quantum algorithm 
which uses k = comp qub_std (e, S) qubits produces at most 2 k different elements A n k (f,j) = 
<fi(j) from G for j = 0, 1, 2, . . . , 2 k — 1 with independent of /. 
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We take an arbitrary quantum algorithm A n ^ with error e qua std (A n) k) < s. For any 
/ G F, we have 

/ 2 *-i \ 1/2 

min ||S(/)-0(j)|| < £p/jW)-M)|| 2 

.=0,1,..^-! ^ y 

Let B(<f)(j),s) = {g E G : \\g — <f>(j)\\ < £ } be the ball in G of center <f)(j) and radius e. 
Obviously diam(B(<j){j),e)) < 2e. Then S(f) E uf^ 1 B ((j>(j) , e) and therefore S(F) C 
U |lo 1 ^(0(j), e). This means that 2 k > n(e, S(F)), and k > log n(e, 5(F)), as claimed. □ 

The essence of Theorem 3.2 is that 5(F) must have a finite e-entropy in order to have S 
solvable in the standard quantum setting. In particular, this means that the closure of S(F) 
must be compact. Otherwise, the e-entropy is infinite and the finite number of qubits is not 
enough to solve S. We summarize this in the following corollary. 

Corollary 3.2. If the closure of S(F) is not compact then S is not solvable in the standard 
quantum setting. In particular, if S(F) is unbounded then S is not solvable in the standard 
quantum setting. 

The unboundedness of S(F) can happen even for problems with relatively small worst 
case complexity as shown in the following example. This example also shows that lower 
bounds based on the e-entropy of S(F) presented in Theorem 3.2 may be quite different 
than lower bounds based on the non-adaptive information complexity in the worst case 
setting presented in Theorem 3.1. 

Example: Unbounded S(F) 

Consider the univariate integration problem for Lipschitz functions, i.e., 

F= {/:[0,1]-R| \f(x)-f(y)\<\x-y\ V x, y e [0, 1] }, 

and S(f) = Jo fit) dt with G = R. 

Since all constant functions belong to F, we have S(F) = R and therefore Ent(e, R) = oo. 
It is well known that the worst case complexity is roughly l/(4e), see e.g., [33], and the linear 
algorithm 

J— * 

with n = \e' 1 /A\ minimizes the worst case error among all algorithms using n function 
values, and has error at most e. Observe that for constant functions, f(t) = c, we have 
= c which may be arbitrary large. 

In the worst case setting with the real number model, the sizes of numbers do not matter 
and do not affect the cost analysis. In the standard quantum setting, the situation is different 
since we can only work on unit vectors and the scaling of numbers does matter. That is why 
we cannot solve unsealed problems in the standard quantum setting. 

In many cases, we may rescale the problem by changing F to a set F such that S(F) 
is bounded and its closure is compact. This idea works for our example as follows. For 
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/ G F, define g(x) = f(x) — /(O) for x G [0, 1]. Then g(0) = and |<7(x)| < x < 1. Hence, 
—a; < g(x) < x and therefore £((?) G [— |, |], and both bounds are sharp. Define 

F = {<? : [0, 1] - R| (?(0) = 0, g G F}. 

Then / G F iff <? G F for <?(x) = f(x) - /(0), and = S(<?) - /(0). 
We now have S(F) = [— \, \\ and therefore 

Ent(£,S(F)) = log e' 1 + 0(1). 

Hence, the number of qubits for approximations of S(g) is now bounded by roughly log e~ l . 
In fact it is sharp, since S(g) can be approximated to within e in the standard quantum 
setting by A n k (q) using roughly n = e -1 / 2 bit queries and k = log e~ l qubits as shown 
in [24]. 

Finally, we may approximate S(f) for / G F by running A n ^{g) for g(x) = f(x) — /(0), 
and computing /(0) + A Hik (q) on a classical computer. Note that the last step on a classical 
computer may involve an arbitrarily large number /(0) which is of no relevance as long as 
we use the real number model of computation. □ 

4 Quantum Setting with Randomized Queries 

Modulo measurements, the standard quantum setting for continuous problems is similar to 
the worst case setting with non-adaptive information. All unitary matrices including queries, 
as well as the number of qubits used by quantum algorithms are deterministic. The potential 
speedup of the standard quantum setting for continuous problems over the worst case setting 
relies on the fact that quantum algorithms with k qubits may use an exponential number 
up to 2 k function values with cost proportional to a small power of k. If 2 k functions values 
are not enough to solve the problem in the worst case setting then the problem remains 
unsolvable also in the standard quantum setting. As we indicated before, this may happen 
even for problems with small randomized complexity. Such examples suggest studying more 
general quantum settings. 

In this section, we describe the quantum setting with randomized queries in which all 
unitary matrices including queries may be randomized. Modulo measurements, the quantum 
setting with randomized queries for continuous problems will be similar to the randomized 
setting with non-adaptive information. We assume that the number of qubits is fixed and 
does not depend on randomization. As we already mentioned in the introduction the exten- 
sion to randomized qubits is left for future study. 

We generalize (6) by allowing unitary matrices Uj as well as the query Qf to be randomly 
chosen similarly as in the randomized classical setting of Section 2.2. That is, we have random 
elements uj distributed accordingly to some probability measure p on Q with p(fl) = 1. We 
stress that p does not depend on / and we will be using the same randomization for all / 
from F. 

First we choose k as the number of qubits, take a random element u, choose a unit vector 
iVfyw) from C 2k as the initial state, and obtain the final k qubit state 

\ipf,u>) = Urtv^QfaUnv-^Qfa ■■■ U^Qf^Uo^llpQ^). (13) 
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For a fixed u, we have the same situation as in the standard quantum setting. That is, 
matrices Uj >u> are arbitrary 2 k x 2 k unitary matrices which are independent of /, and the 
query Qf tU >, which is also a 2 k x 2 k unitary matrix, depends on at most 2 k sample points 
which are independent of / and depend only on u, 

Qf,uj = Qf(ti, u ),f(t2, u ),...,f(t 2 k iU )- 

Hence, Qf,u, is a randomized query depending on at most 2 k function values at randomized 
sample points. In full analogy with the standard quantum setting, the measure p and 
sample points tj jW are the same for all / from F. That is, we use non-adaptive randomized 
information with fixed cardinality at most 2 k . Let 




be the average number of queries used to obtain the final states. 

We then perform a measurement which for a fixed oo is the same as for the standard 
quantum setting. That is, we obtain an index j e {0, 1, . . . , 2 k — 1} with probability P/j> 

depending on the final state \ipf,u>), where Y^jJo~ Pf,j,<*> = 1 f° r a ^ u e ^- As before, the 
dependence on / is only through function values, 



Pf,j,uj - P/(ti, w ),/(fc,„),...,/(t 2fc ,«)j>- 

Knowing the index j, we compute on a classical computer 

for some mapping W : {0, 1, . . . , 2 k — 1} — > G. The algorithm A n>k is called a quantum 
algorithm using randomized queries, or just a quantum algorithm if it is clear from the 
context that we are using randomized queries. 

Analogously to the standard quantum setting, we consider the error of a quantum algo- 
rithm by taking the average performance with respect to both j and iv, see also the appendix 
where the probabilistic error is discussed. That is, the error of an algorithm in the quantum 
setting with randomized queries A n>k is defined by 



1/2 

e q ™(A n , fc ) = sup I / ^p^JSif) - A n M^,j)\\ 2 p(duj)) 
feF \Jn^ Q J 

Observe that if we choose all matrices Uj jU} and Qf,u> as well as n w independently of u, then 
this definition coincides with the error in the standard quantum setting. 

This definition of the error leads to the randomized query complexity defined by 

comp quc - ran (e,5) = min{n : 3A n , k such that e qua - ran (A tjfe ) < e } , 

and to the randomized qubit complexity defined by 

CO mp qub - ran (£, S) = min { k : 3A n , k such that e qua - ran (A Ttife ) < e } . 
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As in the standard quantum setting, we may have a tradeoff between the minimal number 
of queries and qubits. Therefore it would be also reasonable to study the minimization of the 
product or a weighted sum of queries and qubits in the quantum setting with randomized 
queries, however, we do not pursue the issue in this paper. 

It is natural to ask what kind of results can be now achieved and how much we can 
improve the results from the standard quantum setting. We will study these questions in 
the next sections. 



4.1 Lower Bounds on Randomized Qubit Complexity 

For the standard quantum setting, we proved lower bounds on qubit complexity in terms 
of the worst case setting on a classical computer. We now show that lower bounds on the 
randomized qubit complexity can be analogously derived in terms of the randomized setting 
on a classical computer. 

Theorem 4.1. 

comp qub - ran ( £ , S) > log comp inf - ran ( £ ,S). 



Proof. 

We note that any quantum algorithm A n ^ can be regarded as a randomized algorithm whose 
cardinality is at most 2 k . Indeed, for / e F, let pf be a probability measure defined on B x J, 
where B is a measurable subset of Q and J is an arbitrary subset of {0, 1, . . . , 2 k — 1} given 
by 

Pf(BxJ) = / l B (w) J^p /J> p(dw) 
with the characteristic function Ib(u) = 1 for uo E B and 1_b(c<j) = otherwise. Then 

„ 2 fc -l 

eqU a- rEm(Aifc) 2 = gup / J2p f ,,JS(f) ~ A n , k (f,00,j)\\ 2 p(duj) 

feF Jn 

= sup / \\S(f)-A n , k (f,u;,j)\\ 2 p f (d(u;,j)) = e ran (A n , fc ) 2 . 

feF Jnxio,i,...,2 fc -n 



Hence, A njk can be regarded as a randomized algorithm in the randomized classical setting 
whose cardinality is at most 2 k . Since the sample points tj jW are independent of / and 
dependent only on uj, the algorithm A n ^ uses non-adaptive information, and A n ^ applied to 
/ uses randomization with the measure pf. If e qua, ~ va,n {A nk ) < e then e ra,n (A n ^) < s which 
may happen only if the cardinality of A n ^ is at least comp inf ~ ran (e, S). This means that 
2 k > comp inf - ran (e, S), as claimed. □ 

Note that the lower bounds in Theorem 4.1 are not larger than the lower bounds in 
Theorem 3.1 since comp mf ~ ran (£, S) < comp mf_wor (e, S). Furthermore for some problems 
the non-adaptive information complexity in the worst case setting may be infinite whereas 
its randomized counterpart is relatively small. 
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We already mentioned that the integration problem for the class F = F di0 is unsolvable 
in the standard quantum setting and solvable in the randomized classical setting. We now 
provide a proof and find the randomized qubit complexity. The randomized query complexity 
for Fd t Q as well as the integration problem for F^ r for an arbitrary integer r > will be studied 
later. 

Example : Multivariate Integration for r = (continued) 

We show that the randomized qubit complexity for F = F d0 is 

comp qub - ran (e,INT di0 ) = (log e' 1 ) . 

and is achieved by a quantum algorithm which uses of order e~ x queries. 

We know that comp mf_ran (£, INT^o) = @(^~ 2 )- From Theorem 4.1 we conclude that the 
randomized qubit complexity must be at least of order log e v . 

We now provide an upper bound on the randomized qubit complexity. Take the Monte 
Carlo algorithm with m = [4e~ 2 ] , 



1 m 

MC m (i» = - Y,fM 



with uj = [wi,W2,...,w m ] and independent uniformly distributed ojj from [0, l] d . Then 

e ran (MC m )<e/2. 

We now apply the Boolean summation algorithm A* n k of [5] with seven repetitions for 
real functions for which \f{x)\ < 1, see also Section 5. The algorithm A* nk approximates 
MC m (/). It is known, see [11], that the randomized error of this algorithm is bounded by 
C/n, where C is a number independent on /, n and m. Furthermore, the algorithm uses 
k = ©(log m) qubits. We set n — \2C/e] and obtain 



2 fc -l 



^ P/J> [Mc m (/, w )-i; fc (/, w ,j)] 

3=0 



2 < 1 



V/eF. 



Therefore for any / G F we have 



lNT dfi (f)-A* (f, u ,j) 



INT d)0 (/) - MC m (i» + MC m (/,w) - A* k (f,u,j) 



< 



2 k -l 

3=0 

2 k -l 

3=0 

2 k -l |- 

2 ^ P/J> (INT di0 (/)-MC m (/,c)) 2 + {MC m (f^)-K,k(f^^)Y 

3=0 L 



INT di0 (/)-MC m (/,o;) 



< 
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Taking the integral over f2 we conclude that 




Hence, we can solve the integration problem for r = using of order log e 1 qubits and e 1 
randomized queries, as claimed. □ 

Theorem 4.1 states that the number of qubits in the quantum setting with randomized 
queries depends on the non-adaptive information complexity comp mf ~ ran (e, S). Typically, 
comp mf ~ ran (£, S) goes to infinity with e tending to zero, so the number of qubits has to go 
to infinity as well although much slower due to the presence of the logarithm in the bound 
of Theorem 4.1. However, if comp mf_ran (e, S) — oo then the randomized qubit complexity is 
infinity and the problem cannot be solved. We summarize this fact in the following corollary. 

Corollary 4.1. If the non-adaptive information complexity of S in the randomized setting 
is infinity then S cannot be solved in the quantum setting with randomized queries. 

We illustrate Corollary 4.1 by a problem with infinite non-adaptive information complex- 
ity in the quantum setting with randomized queries. 

Example : Multivariate Approximation 

Consider the same class F = F^ r , r > 0, as for the multivariate integration problem. Let 
G = C([0, l] d ) and S = APP d , r : C r ([0, l} d ) -> G be defined by 

APPd, r (/) = /■ 

It is known, see [32] p. 425, that randomization does not help for this problem and that for 
small e we have 

comp inf-ra„ (£)A pp^ o) = comp^"™^, APP40) = OO, 

and for r > 1, 

comp inf - ran (£,APP dir ) = (comp inf - wor (£,APP dir )) = (^ /r ) . 

Hence, for r = we conclude that the approximation problem cannot be solved in the 
quantum setting with randomized queries. □ 

5 Boolean and Real Summation 

Solution of the real summation problem is a basic module used in the solution of many 
continuous problems. The real summation problem is very much related to the Boolean 
summation problem which has been thoroughly studied in the standard quantum setting, 
see [5, 7, 11, 16, 18]. In this section we study the Boolean and real summation problems in 
the quantum setting with randomized queries. 
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5.1 Boolean Summation 

For N a (large) power of two, consider the class of Boolean functions 

F = F N = {f\ /:{0,1,...,7V-1} - {0,1}}, 
and take G — R. The Boolean summation problem is defined as S = B^ '■ F N — > G given by 

1 

3=0 

We want to approximate Bn to within e. Without loss of generality we may assume that 
e > l/(2N). Indeed, if we know A n>k (f) such that \B N (f) - A», fc (/)| < e < 1/( 2A then we 
can recover Bjv(/) exactly since 

« m \NA n ^f) + \-\ - 1 
«iv(/) = • 

This follows from the fact that we know a priori that Z3jv(/) = k/N lor some integer fee [0, iV] 
with being the total number of the true assignments of the Boolean function /. Then 
NA Uik (f) + \ = k + x with 

x = \ + N (A ntk (f) - B N ) e [\ - Ne, \ + Ne] c (0, 1). 

Hence, \k + x] — k + 1, as claimed. 

We first consider the Boolean summation problem in the standard quantum setting. 
The Boolean summation algorithm A* nk of [5] with seven repetitions solves the Boolean 
summation such that 

e^- std (A n , k ) < e, 

using n = Q(e^ 1 ) bit queries and k = 0(log N) qubits. The query bound is order-optimal, 
see [11, 18]. The qubit bound is also order-optimal since it is known, see e.g., [24], that in 
the worst case setting 

comp inf - wor (e, S) = \N(1 - 2e)] G [0, \] 

which is essentially N for small e. From Theorem 3.1 we conclude that the qubit complexity 
is roughly at least log N for small e. We summarize these results in the following theorem. 

Theorem 5.1. The complexities of the Boolean summation problem in the standard quantum 
setting satisfy 

comp^-^Siv) = 6 (e- 1 ) comp qub - std (i37v) = 6 (log N) . 

Furthermore, these bounds are both attained by the Boolean summation algorithm with seven 
repetitions. 

We now consider the Boolean summation problem in the quantum setting with random- 
ized queries. We prove the following theorem. 
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Theorem 5.2. The complexities of the Boolean summation problem in the quantum setting 
with randomized queries satisfy 

comp quc - ran (i3jv) = (s- 1 ) comp qub - r£U1 (£Ar) = 6 (log e' 1 ) . 

Furthermore, these bounds are both attained by the Boolean summation algorithm with seven 
repetitions applied to 



^ in 

with m = 6(e~ 2 ) and with independent uniformly distributed uij from {0, 1, . . . , N — 1}. 



Proo/. 

We first consider lower bounds and start with the randomized query complexity. We use the 
known proof technique of using the average case error as a lower estimate of the randomized 
error. More precisely, take an arbitrary quantum algorithm A n ^ that uses n randomized bit 
queries and k qubits, and consider its randomized error 



e 



2—1 

'""'-'■"(A,,*) 2 = sup /' ^ tfys/) A Hik (f,u;,j)\ 2 p(du;). 
feF N Jn 



We now replace the supremum over / by an average over /. That is, we assume that a 
Boolean function / from Fn occurs with probability pf with non-negative pj such that 
Ylf£F N Pf = 1- Observe that this is a well defined measure since Fn consists of finitely many 
Boolean functions, in fact, we have 2^ functions in F N - Then 



Jn \feF N j=o J 

The following result proved by Papageorgiou in [25] will be needed for our consideration. 
Take the uniform distribution for B N (f), i.e., p f = ( N j )/2 N for all / with B N (f) = j/N. Then 
there are two positive numbers C\ and C2 with the following properties. For any algorithm 
A n ^ in the standard quantum setting with n bit queries, such that n < C\N , and k qubits, 
let pfj denote the probability of obtaining the index j through measurement. Let 

denote the probability of a subset J of {0, 1, . . . , 2 k — 1}. Then it is proved in [25] that 

5>/MO'= \B N (f) - A n:k (fJ)\ > c 2 /n}) > 0.25. 
feF N 
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From Chebyshev's inequality we conclude that 



2 k — 1 

Ep/En/)-^(/-j)I j ^ id)'- 



/gFjv j=0 

We apply the last inequality for an arbitrary algorithm A n ^ from the quantum setting with 
randomized queries and with a fixed u. Here, we use the fact that the algorithm A Ui k(-, 00, •) 
can be regarded as an algorithm from the standard quantum setting. We thus have 

2 k — 1 

E pf E p/^im/) - A>, fe c/>, j)i 2 > 1 ) 2 • 

Since the right-hand side is independent of cu, from (14) we obtain 

Hence, e qua ~ ran (A ntk ) < e implies that n = and 

comp que - rail (B J v) = VI (e- 1 ) . 

To prove a lower bound on the randomized qubit complexity, we use Theorem 4.1. In 
the randomized setting on a classical computer it is known that the randomized complexity 
is of order e~ 2 , see [24]. Then Theorem 4.1 yields 

comp qub - raJ1 ( 1 Biv) = n (e- 1 ) . 

We turn to upper bounds. The idea is the same as for multivariate integration for r = 
which was studied before. That is, we apply the Boolean summation algorithm with seven 
repetitions to the Monte Carlo algorithm mT x X]j=i f( u j) with independently and uniformly 
distributed 00 j over {0, 1, . . . , N — 1}. We stress that this algorithm uses randomized queries 
Qf tUJ and the rest of unitary matrices are deterministic, i.e., Uj tU1 = Uj in (13). 

The same analysis done for r = yields that the randomized error is e. Since this 
algorithm uses of order e -1 randomized queries and log e _1 qubits, we obtain upper bounds 
which match the lower bounds. This completes the proof. □. 

It is interesting to compare the complexities of the Boolean summation problem in the 
quantum settings with deterministic and randomized queries. As we see, the query com- 
plexities are roughly the same in both settings. The qubit complexities, however, are quite 
different. For deterministic queries, the number of qubits depends on the common domain of 
Boolean functions, and we need roughly log N qubits which can be arbitrary large for large 
N. For randomized queries, the number of qubits does not depend on the common domain 
of Boolean functions. It depends on the error parameter through log £ _1 . As we shall see 
in the following sections, there is sometimes an exponential difference between log iV and 
log e -1 . 
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5.2 Real Summation 

We finish this section by a brief note on the real summation problem. We now consider 
/ : {0, 1 . . . , N — 1} — > [0, 1] and want to approximate 

N-l 

SUM Ar (/) = - J^fU)- ( 15 ) 
A known idea is to replace the real number f(j) from [0, 1] by its binary expansion, 

oo 

fU) = 5^2-*/(i,j) with/(t,j) G {0,1}, 

i=i 

We define K = [log e -2 ] and truncate f(j) to if bits. Let 

j=0 i=l 

Clearly, |SUMjv(/) - «S'a'(/)| < 2~ x < e 2 . To obtain a Boolean function we finally define the 
set 

D = {(i,j,p) :i = l,2,...,K,j = 0,l,...,N-l,p = l,2,...,2 K ~ i } 
of cardinality N(2 K — 1) and a Boolean function bf : D — > {0, 1} by 

6(i,j',A;) = 

where is the Kronecker delta. Then 

A - K K 2 K ~ i 

i=l i=l i=l p=l 

Thus 

JV-1 X 2 K ~ i 

skU) = B N2K {b f ) = ^EE^j'-p) 

j=0 1=1 p=l 

is a Boolean summation problem. We compute B N2 K(bf) with error e — e 2 by the Boolean 
summation algorithm with seven repetitions as explained in Theorem 5.2 and obtain A n ^(bf) 
which uses of order e~ l randomized queries and log e _1 qubits. It is easy to check that 
An,k(bf) approximates S(f) with the randomized error at most e. We summarize this in the 
corollary. 

Corollary 5.1. The complexities of the real summation problem (lb) in the quantum setting 
with randomized queries satisfy 

comp quc - ran (SUMjv) = (s- 1 ) comp qub - ran (SUMjv) = (log e' 1 ) . 

Furthermore, these bounds are both attained by the Boolean summation algorithm with seven 
repetitions applied to 

j m 

— y2f(ie,je,Pe) 
l=i 

with m = Q(e~ 2 ) and with independent uniformly distributed (ie,je,Pe) over D. 
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6 Multivariate Integration 



In this section, we consider the multivariate integration problem INT^ ir for the class F^ r 
which we used throughout as an illustrative example. The purpose of this section is to study 
this problem in the quantum settings with deterministic and randomized queries. 

We begin with deterministic queries. For r = the problem is unsolvable. For r > 1, 
sharp bounds on the bit query follow from [24], 

comp quc - std (e, INT dir ) = Q (e~ 1/{1+r/d) ) . 

This bound is achieved by the Boolean summation algorithm with seven repetitions and uses 
of order log e^ 1 qubits, as shown in the previous section. 

Observe that the lower bound on the qubit complexity is of order log e _1 due to Theo- 
rem 4.1 and the fact that the worst case information complexity is of order £~ d/,r . Hence, we 
have 

CO mp qub - ran (£, INT d>r ) = (log a' 1 ) . 

We now turn to randomized queries. The case r = has already been covered and we 
know that we can solve the problem using of order e -1 randomized bit queries and log e~ l 
qubits. 

For r > 1, we use the same quantum algorithm as in [24]. Since the Boolean summation 
algorithm uses the same order of bit queries for the randomized and probabilistic quantum 
errors, we obtain the same upper bounds on the number of bit queries and qubits. 

The lower bound on the number of randomized queries can be derived as in [24] and 
using the results on the Boolean and real summation problems of the previous section. This 
yields that the randomized bit query complexity is of order e~ 1 /( 1 +»7 c The lower bound on 
the number of qubits follows from Theorem 4.1 and the fact that the randomized complexity 
is of order e -2/(i+2r/«0 . 

We summarize these results in the following theorem. 

Theorem 6.1. Consider the multivariate integration problem INT djr for the class F^ r . 

• Letr = 0. 

— In the quantum setting with deterministic queries, we have 

comp quc - std (£,INT di0 ) = oo, 
comp qub - std (e,INT, i0 ) = oo. 

— In the quantum setting with randomized queries, we have 

comp^-^MNT^o) = e^- 1 ), 
comp qub - ran ( £ ,INT d) o) = G^oge- 1 ). 

• Let r > 1. 
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— In the quantum setting with deterministic queries, we have 

comp quc - std (e,INTV) = Q (e~ 1/{1+r/d) ) , 
comp qub - std (e, INT d , r ) = (log e' 1 ) . 

— In the quantum setting with randomized queries, we have 

comp qu ™(£,INT dir ) = ( e -V(i-H-A0) ) 
comp qub - ran (e,INT d , r ) = (log e' 1 ) . 

Hence, for r = we see a big difference between the two settings for the multivariate 
integration problem, whereas for r > 1, the two settings lead to the same order of bit and 
qubit complexities. 

7 Path Integration 

Path integration can be regarded as integration of functions of infinitely many variables 
or, more formally, as integration over some class of functions; for more information and 
references see [35]. Path integrals occur in quantum physics, chemistry and mathematical 
finance. They are also the solutions of certain differential equations and mathematical finance 
problems. 

Here we consider a specific example of path integration studied in [35]. We take the space 
X := C([0, 1]) of continuous functions defined on [0, 1] with the norm ||a;|| = max tG [ 0i i] \x(t)\. 
The space X is equipped with the classical Wiener measure w for which 

/ x(t)w(dt) = Vie [0,1] and / x(t)x( u)w(dx) = min(t,-u) Wt,uE [0,1]. 
Jx J X 

We consider the class F of real valued w-integrable functions / : X — > M which are bounded 
and satisfy a Lipschitz condition. More precisely, let the norm of / be given by ||/|| = 
snp xe x \f( x )\- Then the class F is defined as 

F = {f : H/ll < 1, \f(x)-f(y)\ < \\x-y\\ L2m]) Vx,yeX}. 

Let G = R. The path integration S := PATH is given by 

PATH(/) = / f(x)w(dx). 
Jx 

We first consider the standard quantum setting. It was shown in [35] that we can compute 
an e- approximation for path integrals from the class F with probability | using of order e~ l 
bit queries and e~ 2 log e -1 qubits. The bound on the number of bit queries is sharp in the 
sense that for any positive a it cannot be smaller that e l ~ a as e goes to zero. The sharpness 
of the number of qubits was not discussed. These bounds are obtained by reducing the path 
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integration problem to the summation problem which was solved by the Boolean summation 
algorithm. 

We consider the error defined by (9) in the standard quantum setting. Due to the fact 
that the Boolean summation algorithm enjoys optimality properties also in this setting, we 
conclude that the same bounds as above also hold for the error (9). Furthermore, the bound 
on the number of qubits is sharp since the worst case information complexity is of order 
£ -ce 2 £ or some positive c as proved in [6]. Then Theorem 3.1 yields that the number of 
qubits must be of order e~ 2 log e^ 1 . 

We now consider the quantum setting with randomized queries. It was proven in [37] 
that the non-adaptive information complexity in the randomized classical setting is of the 
form 

comp inf - ran (£,PATH) = ( e -2(i+°(i))) as £ _, . 

This and Theorem 4.1 yields that the number of qubits is at least of order log e~ x . 

We now show that in the quantum setting with randomized queries, we can solve the path 
integration problem by using of order e^ 1 bit queries and log e" 1 qubits. The space X can 
be embedded in the Hilbert space L 2 ([0, 1]) for which the embedding Im : X — > L2QO, 1]), 
Imx = x for all x G X, is a continuous linear operator. Let v = wlrn -1 be a zero mean 
Gaussian measure on L2QO, 1]). Then the covariance operator C v of the measure v has 
eigenpairs, C v r\i = A^, where 

r)i(x) = V2 sin - - irx^J , A* 

As in [35] we first approximate PATH(/) by 

INT d (/) = / f d (t)fjL d (dt), 

where 

f d (t) = f {Im' 1 (hrn + t 2 7] 2 + ■■■ + t d r] d )) 

and \x d is the d dimensional Gaussian measure with the mean zero and variances A«. That 
is, its density function is of the form 

exp(-t 2 /(2A 1 ) tl/{2X d )). 



7T 2 (2i - l) 2 ' 



(27r) d /2VA 1 A 2 ---A, 

In [35], it is proved that for d = Q(e- 2 ) we have |PATH(/) - INT d (/)| < e/3 for all f e F. 
The integral INT^(/) can be approximated by the Monte Carlo 



with iid points tj distributed according to the measure /j, d . Note that for / e F, we have 
\fd{tj)\ < 1 and clearly the variance of f d is bounded by 1. Therefore for n = \9e~ 2 ], the 
randomized error of approximating INT d (f) is at most e/3, and the randomized error of 
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approximating INT(/) is at most 2e/3. Finally, it is enough to use the Boolean summation 
algorithm to approximate the last sum with the randomized error e/3 which can be done 
with of order e' 1 bit queries and log n = 0(log qubits. The randomized error of 
approximation INT(/) is at most claimed. 

This and the previous lower bound on the number of queries yield that the randomized 
qubit complexity of path integration is of order log e~ x . For the randomized query complexity 
we have so far an upper bound of order e~ l . We can get a lower bound by applying the same 
proof technique as in Theorem 3 of [35]. That is, the path integration problem is reduced 
to the real summation problem for which we use a lower bound presented in Corollary 5.1. 
This yields that 

lim£ 1 - Q comp que - ran (£,PATH) = oo Va G (0,1). 

We summarize these results in the following theorem. 

Theorem 7.1. Consider path integration equipped with the Wiener measure for the class F 
of Lipschitz functions. 

• In the quantum setting with deterministic queries, we have 

comp quc - std (£,PATH) = 0(£- 1+o(1) ), 
comp qub - std (e, PATH) = (e' 2 log e' 1 ) . 

• In the quantum setting with randomized queries, we have 

comp quc - ran (£,PATH) = e(e _1+o(1) ), 
comp qub - ran (e, PATH) = (log e' 1 ) . 



The essence of this theorem is that for path integration we have an exponential improve- 
ment in the number of qubits in the quantum setting with randomized queries whereas the 
number of queries remains roughly the same in both settings. We stress that the optimal 
bounds for bit queries and qubits are both attained by the same quantum algorithm based 
on the Boolean summation algorithm. 



8 Appendix: Probabilistic Errors 

We briefly indicate what kind of results are possible if one studies probabilistic errors in the 
randomized classical setting and in the quantum settings with deterministic and randomized 
queries. 

We begin with the randomized classical setting. Instead of the randomized error (3) we 
now consider the probabilistic error of the algorithm A n which is defined by the worst case 
performance with respect to / and the worst case performance with respect to uo modulo a 
set of measure 5 for some (usually small) 5 G (0, 1). That is, 

e^(A n ,5) = sup inf sup \\S(f) - A n (f,u)\\. (16) 

fzF Bcn,p(B)<s we n\s 
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From Chebyshev's inequality we have 



e^\A n ,5) < ^e ran (^n). (17) 

Better estimates with respect to 5 are available under additional assumptions on S. In any 
case, the dependence on 5 is quite modest and everything depends on the randomized error 
e ran (A n ). This is probably why the probabilistic error (16) has not been as widely studied 
as the randomized error (3) for continuous problems on a classical computer. 

The probabilistic error yields the (non- adaptive) information complexity defined by 

comp inf - ran (e, 5, S) = min { n : 3 A™ A such that e ran (A° ad , 5) < e} . 

Clearly, (17) implies that 

comp inf - ran ( £ , 5, S) < comp inl - ran {eVS, S), 



We now turn to the standard quantum setting. The usual way of defining the error in 
this setting is analogous to the probabilistic error. That is, the probabilistic error of A n>k is 
defined as the smallest a for which 

\\S(f)-A n , k (f,j)\\ < a 

holds with probability at most 1 — 5 with respect to j for every / from F. Here, S G (0, 1). 
This definition can be formalized as follows. For / e F and an arbitrary subset J of 
{0, 1, . . . , 2 k — 1}, let f-if(J) = J2jejPfd De a measure of J. Then the probabilistic error of 
A nt k in the standard quantum setting is 

e^- std (A n , k ,5) = sup min max - A n , k (f, j)l 

feF J:M/(J)<5 ie{0,l,-,2 fe -l}\J 

For some operators S, such as linear functionals, it is typical to take, say, S — |, and 
obtain a quantum algorithm working with probability 1 — 5 by repeating several times the 
quantum algorithm working with 5 — | and by taking the median as the final result. If the 
number of repetitions is large enough we can boost probability of success to 1 — 5. Details 
can be found in [7]. 

The probabilistic query complexity in the standard quantum setting is defined as 

comp quc - std (e, 5, S) = min { n : 3 A n , k such that e qua - std (A njfc , 5) < e } , (18) 
and the probabilistic qubit complexity in the standard quantum setting is defined as 

comp qub - std ( £ , 5, S) = min { k : 3A n>k such that e qua - std (A„ ife , 5) < e } . (19) 

Example : Multivariate Integration (continued) 

Assume that r > 1. For F = Fd jr , it has been proven by Novak, see [24], that for 8 = \ the 
minimal error of quantum algorithms A n>k is of order n~ 1 ~ r ^ d , and is achieved by an algorithm 
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that uses of order log e~ x qubits. The idea of the proof was to reduce the integration problem 
to real and then Boolean summation and apply the the Boolean summation algorithm of 
[5]. This algorithm for the summation of N terms with n queries, n <^ N, has probabilistic 
error of order n" 1 which is optimal due to [18]. 

This implies that the query complexity comp quc ~ std (e, |, INT djf ,) is of order £~ 1 /( 1+r / d ). 
For d much larger than r, we thus obtain roughly a quadratic speedup over the randomized 
setting, and an exponential speedup over the worst case setting. 

For arbitrary 5, we can use roughly log 5" 1 repetitions of the algorithm used for 5 — | 
and take the median of computed results as the final result, again see [7, 21]. This implies 
that comp quc - std (£, 5, INT d>r ) is of order e-Va-H-AO i og §-\ 

For r = 0, it is easy to see by the same argument which we used for the randomized errors, 
that the integration problem cannot be solved for the probabilistic error in the standard 
quantum setting. □ 

As before, Chebyshev's inequality yields 

comp quc - std (e, 5, S) < <T 1/2 comp quc ~ std (e, S) 
comp qub - std (e,5,5) < (T 1/2 comp qub - std (£,S). 

Again the dependence on 5 can be improved for some S. Note, however, that even for 
general S, the dependence on 8 is quite weak. 

We now show lower bounds on the probabilistic qubit complexity in terms of the non- 
adaptive information complexities in the worst case and randomized settings as well as in 
terms of the ^-entropy. 

Theorem 8.1. 

comp qub - std (e,(5,5) > log comp inf - wor (2e, S) V£e(0,±), 

comp qub - std (e,(5,5) > log comp inf - ran (£, 5, S), 

comp qub - std ( £ , 5, S) > Ent{e, S{F)). 



Proof. 

(1) To prove the first inequality, we take a quantum algorithm A n ^ which uses the minimal 
number of qubits k = comp qub - std (e, 5, S) with e qua - std (A„ ife , 5) < e. We have 

£ > e qua-std (A 5) = gup m j \\ S (f)-A n , k (fj)\\ (20) 

feF J:n f (J)<8 jG{0,l,...,2 fc -l}\J 

= sup max \\S(f) - A n<k (f,j)\\, 

/GFje{0,l,-,2 fc -l}\J(/) 

where J(f) is a subset of {0,1,..., 2 k — 1}, < 8, for which the corresponding 

minimum is attained. Such a set exists since we have finitely many such subsets, however, 
J(f) is not necessarily unique. Let 

M(f) = {0,l,...,2 fc -l} \ J(f). 

Clearly, n f {M{f)) >1-S. 
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For an arbitrary / G F, we take two functions fi and f 2 such that N(fi) = N(f 2 ) = N(f). 
Note that the measures /i^ and jif 2 are the same. For 5 < |, there exists an index j* which 
belongs to M(/i) fl M(/ 2 ). Indeed, otherwise the sets M(/i) and M(/ 2 ) would be disjoint 
and 

1 > /x/(M(/i) UM(/ 2 )) = / x / (M(/ 1 ))+///(M(/ 2 )) 

= iM h (M(fi))+^(M(f 2 )) > 2(1-5) >1 

is a contradiction. For this index j*, we have a = A n ^ k (fi, j*) = A U:k (f 2 ,j*). From (20) we 

^ e>5 (WO " «ll + Wa) - a||) > |||5(A) - 5(/ 2 )||. 

Hence, 

sup sup ||5(/i)-5(/ 2 )|| < 2e, 

/€F fi,f2eF,N(f 1 )=N(f 2 )=N(f) 

and the rest of the proof is the same as in the proof of Theorem 3.1. 

(2) To prove the second inequality, we compare the qubit complexity to the randomized 
non-adaptive information complexity. Observe that any quantum algorithm A n>k may be 
regarded as a randomized algorithm which uses non-adaptive deterministic information of 
cardinality at most 2 k with randomized elements uo G {0, 1, . . . , 2 k — 1} taking values j with 
probability pfj. Furthermore, the probabilistic error of A n ^ is exactly the same as the 
error in the probabilistic randomized setting. Therefore, e van (A n) k, 5) < e can hold only if 
the cardinality of A n ^ is comp inf ~ ran (£:, 5, S). This means that 2 k > comp inf ~ ran (e, 5, S), as 
claimed. 

(3) To prove the third inequality, observe that we now have for any / e F, 

.max } \\S(f)-cf>(j)\\ < e, 

where the subset M(f) of {0, 1, . . . , 2 k — 1} is defined as above. Hence, 

min \\S{f) - 4>{j)\\ < e, 

j=0,l,...,2 fc -l 

and the rest is as in the proof of Theorem 3.2. □ 

We now turn to the quantum setting with randomized queries. Instead of (9), we consider 
the probabilistic error of A n>k which is defined as the smallest a for which 

\\S(f)-A n , k (f,ujJ)\\ < a 

holds with probability at least 1 — 5 with respect to j and uj for all / from F. More precisely, 
as before, for J e {0, 1, . . . , 2 k — 1} we define the measure of J by /i/ jW (J) = ~^2j e jPf,j,uj- 
Then the probabilistic error of A n ^ k in the quantum setting with randomized queries is 

e qua-ran^ ^ = j nf m&x \\S {f) - A Uik (f , UJ , 

feF BeCl,p(B)<8 uen \ B J:n flU (J)<6 je{0,l,-,2 fc -l}\J 
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Observe that if we choose A n ^ independently of to then, modulo measurement, everything 
will be deterministic and the last definition coincides with the probabilistic error in the 
standard quantum setting. 

The probabilistic query/qubit complexity in the quantum setting with randomized queries 
are defined, analogously as in the standard quantum setting, by minimizing the number of 
queries/qubits needed to find a quantum algorithm whose probabilistic error is at most e. 
That is, the probabilistic query complexity in the quantum setting with randomized queries is 
defined by 

comp quc - ran (e,5,S) = min { n : 3A n , k such that e qua ' ian {A ntk ,5) < e } , 

and the probabilistic qubit complexity in the quantum setting with randomized queries is 
defined by 

comp quh ' ian (e,S,5) = min { k : 3A n<k such that e qua - ran (A„ ife , 5) < e } . 
We now show that Chebyshev's inequality implies 

e q "\A n , k ,5) < 5' 1 e q ™(A n , fc ). (21) 

Indeed, let 

2 fc -l 

E(f,u>) = J2Pf^\\S(f)-A n , k (f^J)\\ 2 - 

j=0 

Hence, 

e*"\A ntk f = sup / E{f,u)p{dw). 
feF Jn 

Define the sets B(f) and J(f,Ld) by 

n \ B(f) = {lu : E(f, oo) < 8~ ^-^(A^) 2 } , 
{0,l,...,2 fc --l}\J(/,c) = {j: \\S(f) - A n , k \\ 2 < S^EiM}. 

Chebyshev's inequality tells us that p(B(f)) < 5 and fif tLU (J(f,oj)) < 5. Then 

e^™(A tjfc ,5) 2 < sup sup 6- 1 E(f,u) < r 2 e q ™(A lifc ) 2 , 
feF wen\B(f) 

as claimed. 

The probabilistic and randomized query and qubit complexities are related. From (21) 
we have 

comp quc - ran (£,5,S) < comp qu ™ (e5,S), 
CO mp qub - ran (£, 5, S) < comp (iuh - ran {e5,S). 

It is also easy to see to check that 

comp qub - ran (£, S) > log comp inf - ran (e,.S), 
comp qub - ran (e, 5, S) > log comp inf - ran (5, 5, S). 
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